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1.0 INTRODUCTION 

Over the past twelve years since 1964, 1 have been called upon to perform tests 
to determine the radar (lidar) cross section of several satellite retrodirector 
arrays used in conjunction with laser ranging experiments performed by NASA. 
These satellites have included BE-B, BE-C, GEOS I{A), GEOS n(B), GEOS IH(C) 
and TIMATION m. As a result of this testing I have become familiar with many 
of the problems of designing and testing cube corner arrays. Much of the theory 
behind these tests has not previously been published. I have therefore collected 
many of my notes and assembled them into this document for the reader interested 
in cube corner array performance and testing. It is hoped that this paper will 
explain many of the questions left unanswered in a previous paper, ’’ Design of 
Retrodirector Arrays for Laser Ranging of Satellites ," X-723-74-122. This 
previous paper was Intended as an introduction to cube corner array design 
and discussed arrays made up of diffraction limited cube corners. In practice 
the diffraction limit is seldom obtained due to manufactuidng problems. There- 
fore in dealing with most arrays we have a more difficiilt problem in determining 
areas and gain functions . 

The paper starts with section describing the satellite coordinate system 
necessary to describe the location and orientation of each cube comer in the 
array. The second section describes the method of optical testing. This is 
followed by sections defining the gain function, computational methods for 
deriving the gain function and experimental values for it. The velocity aberration 
is derived as a function of satellite orbit, a complete method for cross section 
evaluation is described, and finally the radar equation is described. 

The methods used in this paper are analytical and do not employ computer 
techniques. This approach has the advatage of giving good insight into the 
problem. However, a computerized analysis is required to perform complete 
analysis. Several computer methods have been utilized by various authors 
with varying degrees of success, but, none of their programs have been documented 
nor are any of the methods capable of accounting for manufacturing errors. 

The author is presently working on a computerized analysis technique for 
cube arrays which will be partly based on previoxis programs, but will extend 
the analysis to include all of the effects described in this paper. 
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2.0 SATELLITE COORDINATE SYSTEM 

Most of the satellites used for laser ranging targets to date have been either 
gravity gradient or otherwise stabilized so that one axis of the satellite remains 
constantly pointed at the earths center. For a laser tracking mission it is de- 
sirable to have the lidar cross section independent of the azimuth angle of the 
satellite viewed from the tracking station. Therefore the retro-director array 
is normally made as nearly as possible, a surface of revolution with constant 
cross section for all azimuth aspects. Before continuing in the diseussicai of 
retrodirector arrays, we must set up a definition of our frame of reference. 

The position and orientation of cube corners in an array are normally specified 
in a spherical coordinate system with a polar axis defined by the line between 
the center of gravity of the spacecraft and the center of the earth as shown in 
Figure 1. The zero reference for the polar angle 9 is the nadir of the satellite, 
and azimuth angles when required are defined relative to the velocity vector of 
the satellite with azimuth angles increasing in a counter-clockwise direction 
when viewed from earth. If the satellite spins about the polar axis, an arbitrary 
azimuth direction related to the structure of the satellite must be selected and 
a time relationship between this azimuth and ihe velocity vector established. 
However, in the case of spinning satellites, the designer normally msdces the 
array rotationally symmetric so that the cross section becomes a function of 
polar angle 6 alone. In the most general case, a cube corner location may be 
specified by three position coordinates (R^ » orientation coor- 

dinates ( 0 ^, 4 >^ ) which specify the direction of the outward normal from the prism 
entrance pupil. The direction of the incident laser radiation may be specified 
by two angles with respect to the satellite coordinate system { 6 ^, 4 )^). 


3.0 OPTICAL TESTING 

In testing a satellite retrodirector array, we must determine the effective 
reflective area, and the gain as a function of the incident laser radiation 
coordinates. This is normally done as shown in Figure 2. The array is placed 
in the beam of an autocollimator and the total energy reflected from the 
array is compared to the energy reflected from a flat of known area and 
reflectivity. From this measurement we can calculate the effective area of 
the array. It should be pointed out that the optical losses in the cube corners 
are taken into accoimt in this measurement method, and therefore no further 
efficiency corrections are necessary. The second part of our test is to scan 
the retrodirected intensity pattern in the collimator focal plane to determine 
the far field point spread function. This is normally done vath a pinhole on an 
XY coordinate stage mounted in the collimator focal plane. Radiation passing 
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through the pinhole is detected by a photomultiplier. Both of these tests must 
obviously be done for all possible values of 6^ and in the case of arrays which 
are not azimuthally symmetric for all Examples of the results of the re- 
flective area and for field pattern tests are shown in Figures 3 and 4. 


4.0 POINT SPREAD FUNCTION 

We can compute the lidar cross section in the following manner. The far field 
pattern is examined and an equation for intensity derived for it (Fig. 5) 


1(0, 7?) = HA^P(0, 7?) 


( 1 ) 


where H is the incident irradiance, A ^ is the effective area and P(0 , ) is the 
point spread function which describes the relative intensities in the far field. 
Since energy must be conserved 


HA 

e 




1(0, 7?) 0d0d?7 


= HA 


e 



P(0, 7]) 0d0d7] 


(2) 


Therefore the definite integral must equal unity. The definition of cross section 
is 


Cr(0, Tj) = 4 tT 

Substituting equation 1 into 3 we obtain 



c^(0, v) = 4t7 


HA^P(0,' 7 )) 
H 


= 4ttA^P(0, tj) = G(0, 7j) Ag 


(3) 


(4) 


Note that H, the irradiance has dropped out of the equation so that in our measure- 
ments it is not necessary to maintain a calibration of the irradiance level. 
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In making the photomultiplier scans of the far field patterns, there are many 
unkn own experimental parameters which make it difficult to measure absolute 
values of intensity. Therefore relative intensities are plotted to an arbitrary 
ordinate determined by the experimental parameters. The intensities are then 
normalized by dividing their values by the ordinate at the (0, 0) position of the 
far field. From these curves, a suitable point spread function P(>/j, 17) must be 
inferred. 


5.0 COMPUTATIONAL METHODS FOR DERIVING THE POINT SPREAD 
FUNCTION 

The most obvious method for deriving the point spread function P(0, rj) would be 
to obtain the Fraunhofer diffraction pattern from the two dimensional Fourier 
transforms of the cube corners in the array. The point spread function could 
then be obtained by summing the amplitudes, squaring the sum and taking the 
absolute value. The result would be proportional to intensity and we could obtain 
P(0 , T 9 ) by a suitable normalization constant. 

This has been done for an equilateral hexagon aperture by Smith and Marsh* vflio 
obtain the following relation 


I 

1 

s in /x_ 
s i n 7^ - 

s i n 

+ 5 i n 7^ 








(5) 


where 

k - 2v/ \ 

X = kl 
Y = km 
= X/3 ± 

= X ± 2Y/-/J 

and ■i and m are the tangents of the rectilinear angles in the far field which for 
small angles are equal to angles themselves. 


* 


R.C. Smith and James C- Marsh, Journ. Opt. Soc. Am. 64, p. 798, June 1974. 
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The comparison of this equation against the measured intensity profiles shown 
in Figure 4 clearly indicates that this does not provide a satisfactory explanation 
of the actual intensity profiles. 


In order to explain the actual patterns we must take the random manufacturing 
errors into account. From " The Scattering of Electromagnetic Waves from 
Rough Surfaces " by P. Beckmann and A. Spizzichino (Mac Millan 1963) we can 
treat the cube corner as a rough surface at normal incidence with the following 
equation 


Kpp*y 


e"g 



exp 



( 6 ) 


where 


> 


g 


Ratio of mean scattered power in the direction specified 
divided by the power in the specular direction for a perfectly 
flat surface i.e. I(i/», t])/\q (p. 29)* 

= cr^ for normal incidence (p. 82)* 



= Fraunhofer diffraction term (p. 82)* 

T = Correlation distance <<-[>2 « ^ 81)* 

F = 1 = Geometry factor at normal incidence (p. 25-27)* 
A = Surface Area (p. 78)* 
small angles 

Therefore the equation reduces to for small angles 


!('/'. ri) _ 

Iq 


pI(.P> i) + 


e"(k'AT/2)^ 

A 


( 6 ) 


Note that the subscript o refers to the quantities for a perfectly smooth surface. 


Refers lo page numbers where definitions in Beckman and Spiz2ichino are found. 
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Some interesting results are that two terms, one for diffraction, and one for 
scattering occur. While the diffraction term is a function of 4^ and 17, the 2^2 
scattering term depends only on . The intensity is reduced by a factor e 
by the randomness of the surface; however this is partially compensated for 
by the scatter term. Since « A the scatter term is much smaller than the 
diffraction term at V' = 0 but will predominate as 4^ grows and the diffraction 
term diminishes. Note also that the scatter term grows as A decreases which 
is the same effect as seen for diffraction. Therefore as the cube corner is 
tilted we would expect a non- symmetrical spreading due to diffraction, plus a 
symmetrical spreading due to increased scatter. 

The random errors in the dihedral angles cannot be taken into account by the 
above equation since they violate the condition that T^ « A. However for an 
array of many cubes with random errors we would expect a Gaussian far field 
pattern if diffraction and scatter could be neglected. 


6.0 EXPERIMENTAL POINT SPREAD FUNCTIONS 

To this point our discussion has not led to a definite conclusion as to what point 
spread function we should expect. There is a strong tendency to expect a dif- 
fraction pattern, or a normal distribution due to random errors depending upon 
what magnitude of fabrication errors are assumed. However, tests have not 
shown this to be true. Rather an equation of the form 


l(s^. v) _ g-cv-zr (7) 

1 ( 0 , 0 ) 


where 


and 




b2 

1 - cos^ 7j 



appears to provide the best fit to the data.* 


For planar arrays 


6 



This equation predicts a point spread function whose iso-intensity contours are 
ellipses with eccentricity e . 

If we let 


P(^, ,) = 

iTih^ 

Then it may be shown that 

Jo Jq 2vb‘^ 

\ 

as required by conservation of energy. The ellipticlty of the pattern is caused 
by the angle of incidence of the radiation at the cube corner. By our definition 
of terms the radiation is incident in the plane containing v = tt/2 so that the long 
dimension of the intensity pattern will always be in the 17 = tt/2 or 3 tt/2 direction 
with 


( 8 ) 


( 9 ) 


( 10 ) 

where is the incidence angle. Substituting € =1.16 \/^ gives a good fit over 
the range ( 9 , = 0 to =77/4 which includes most angles of interest. Therefore 


/I - 1.35 0. 

V. &i) = i. 

27Tb2 


( 11 ) 


where 


1 - 1.35^^ cos^ 17 

The value of c/b can be foimd by examining the point spread functions at normal 
incidence (6. =0). Since only the ratio c/b enters the equation it is not necessary 
to find them independently simplifying. 


77 , 0 .) = 


1/1 - 1.35 0. 


, -/1-1.35 0. 
p-'e ‘ 


2 a 

COS & 

n 


p'l> 


where p = c/b 


2tt 


( 12 ) 



It should be noted that if the cube corner is not rotationally symmetric, the 
point spread function will depend upon the azimuth of incidence <?!>. as well as 
. For the Hexagonal cube corners used on most satellites the minor vari- 
ations due to are neglected. 

Values of p for several satellite arrays determined by analysis of the far field 
intensity patterns are shown in Table 1. 


7.0 THE VELOCITY ABERRATION 

Up to this point tiie need for a two dimensional point spread function has not been 
made clear. The need arises because of the velocity aberration effect which is 
a function of orbital position and the direction of the velocity vector with respect 
to the transmitter/receiver axis. This effect causes the receiver to be located 
at various points in the far field pattern as a fiinction of the orbit and velocity 
vector parameters. 

If a satellite equipped with a cube corner array is moving with a velocity V' 
normal to the transmitter/ satellite axis, a Bradley or velocity aberration effect 
causes the reflected beam pattern to be angularly displaced by an amount 



(14) 


where c is the speed of li^t. Therefore, the reflected radiation pattern is not 
centered on the transmitter, but is displaced by Hie angle xp in the plane con- 
taining tlie velocity vector and the transmitter/satellite axis. 

However, for a satellite in circular orbit, tiie velocity vector will only be normal 
to the transmitter /satellite axis when the satellite is at zenith. The more gen- 
eral case is shown in Figure 6. This figure shows the satellite at zenith angle 6>, 
and altitude h. The satellite and transmitter are in the YZ plane with the trans- 
mitter located at the point (o, o, R^) where R^ is the radius of the earth (6370 Km). 
For purposes of our calculation the earth will be assumed spherical and the orbit 
circular. The satellite is located at the position (o, msin <9 , r cos ^ ). Under 
these conditions the slant range p is 


p - rv - cos 9 
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where 


V = Vl - (R"/r)^ sin^ d 
r = + h 

h = Satellite altitude 

6 - Angle between z axis and slant range p 

Ilie angle between the slant range p and the radius r is denoted as y and is 
equal to 


7=_lsin0 (16) 

r 

The orientation of the satellite reference system to the satellite/transmitter 
axis can be specified by a polar coordinate system centered around the p axis 
with the angular distance from the axis and v its azimuth when viewed in 
the satellite transmitter direction. The zero reference for t) win be selected 
as the -X direction. Therefore a gravity gradient satellite will have its gravity 
gradient axis oriented in the direction 

^ = 7 

Tj - Ztj/'I 

Azimuth angles are measured counterclockwise. 

hi order to define the velocity vector, we must establish its direction and mag- 
nitude. Since the vector is always in the plane normal to the radius r we need 
only specify its direction by an azimuth which we will call oi again referenced 
to the -X direction. 

The velocity normal to the p axis is 


V' 


= V /cos^ CO + s in^ cj 


( 16 ) 


and its direction is 


Tj = tan~^ [v tan oi] 


( 17 ) 
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Two special cases are of interest. The first occurs when w = o or vr. Under 
this condition 


V' = V and 77 = 0 ( 18 ) 

This condition occurs at the highest point in a satellite pass and corresponds to 
the maximum velocity aberration. Note that maximum velocity aberration can 
occur at any zenith angle and that a pass through zenith is not required. The 
second case of interest occurs when co = 77/2 or 3vr/2. When this occurs 


V' = vV 


For satellites in high orbits v « 1 and V’ » V. 

This indicates that for satellites in high orbits, the velocity aberration is nearly 
constant and independent of both co and 0 . 

In contrast, for low orbiting satellites, the velocity V' can vary from V at zenith 
to nearly zero at the horizon. 

The absolute value of velocity aberration is shown in Figure 7 as a function of 
satellite altitude. Data for this curve are shown in Table 2. The variation in 
velocity aberration with d and w is shown in Figure 8 for a 1000km orbit, 
while Table 3 provides values for several orbit heights. 


8.0 EVALUATING CROSS SECTION 

We now have all the equations necessary to calculate the effective cross section 
of a retrodirector array in orbit. The procedure is as follows. Using the data 
on effective area and point spread function, compute cross section using Equa- 
tion 4. The result will be an expression involving four parameters. 


cr = f ($j^, <pj^, Tj) ' ^ 

Next, compute <P^ from Equation 16 (< 9 ^ = 7 ). If the cross section is a function 
of compute 0^. For non-rotating satellites 4>^ = co - it/2. For rotating 
satellites a time fimction must be computed for relating to . Using Equa- 
tions 14, 16, and 17 compute V' > and 77 . Substitute these vali’es of the four pa- 
rameters into Equation 20 to obtain the cross section. Note, that in order to 
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solve for a- we must have the altitude of the satellite, the zenith angle 6 and 
the velocity vector azimuth o> as given data. 

We now have obtained the cross section of the satellite as a function of its orbital 
position. However in order to use cross section we must apply the radar equa- 
tion. Values of cross section are shown for most of the arrays in orbit in Table 5. 


9.0 THE RADAR EQUATION 
The classical radar equation is 


EtGtGi^W 

(4tt)3 r-» 


where 


( 21 ) 


S = Received signal Energy 
E^ = Energy generated by transmitter 
= Antenna gain of transmitter optics 
Gjj = Antenna gain of receiver optics 
A. = Wavelength 

o- = Radar (lidar) cross section of target 
R = Range of target* 

In practical operational conditions a loss factor for the transmission of the 
receiver and transmitter optics and a factor for atmospheric trans- 
mission must be included. Equation 21 becomes 


(4tt)3 R4 


( 22 ) 


Of the parameters in Eq. 22 some are dependent upon the lidar system used, 
some depend upon the propagation path, and one depends upon the target. In 
practice, a critical value of S dependent upon the receiver system and noise 
sources determines the required value of the right hand side of equation 22. 


(47t)^ R“* 

p was used for target range in Section 7. 


( 23 ) 
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If we define three parameters as follows 


P = E_G_G„Trt^^ meter = Stat ion Parameter 

s T T K V 


P = 1 meters'"'^ - Path Parameter 

^ (4t7)3 r4 


P^ ^ o- meters^ = Target Parameter 


then Eq. 23 simplifies to 


P^PpPp ^ 1 (24) 

The reason for defining the parameters in this manner is that each of these 
three parameters is normally under the control of a different group. The track- 
tng station operators can control but is usually controlled by the orbital 
requirements of the spacecraft. The final parameter P^ or simply a- must be 
optimized by the cube comer array designer to make the triple product as large 
as possible but in any case greater than unity. 

We will now analyze the three parameters. 

The required signal S is determined by the number of photoelectrons detectable 
by the detector circuitry. This number N is presently set at about 100 but can 
vary depending upon ranging accuracy requirements and detector noise. Signal 
energy is related to N by the followii^ relation 

S _ N (25) 


where 

hv = Energy of a photon 
7] = Detector quantum efficiency 
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For the most common lasers (Ruby and Doubled Nd - Yag) 

= 5x 10‘i« Joules (\ = 0.69 x 10‘® m) 

(3.8 X 10"^®) 100 

Sj, = ~0~20 “ ^ lO’i® Joules (k = 0.53 x 10“® m) 

The antenna gain is controlled by the 1/e ^ power beamwidth desired (Assuming 
a Gaussian beam profile) by the following equation 


G.J. = 32 / 


(26) 


Receiver gain may be obtained from the classical formula (D^ = receiver 
diameter) '' 


G 


R ~ 



(27) 


System losses must be obtained by a careful analysis of the optical losses 
in the transmitter and receiver optics. 

For the laser tracking stations presently in use by NASA/GSFC the value of the 
station parameter is 


= 5.11 X 10^2 meters^ = 227.1 dB meters^ (28) 

where the following station parameters were used 

E.J. =1.0 Joules 
0^ = 5 X 10~‘* radians 

Dj. = 0.51 meters 
k = 0.6943 X 10“® meters 
Tq = 0.078 
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The path parameter depends upon two variables, satellite altitude and zenith 
angle. For a circular orbit of altitude h 


R =: /(R^ + h)^ - (Rg sin cos 6 (29) 

where = Earths Radius = 6.37 x 10® m. The atmospheric transmission may 
be computed by the following relation 


T 

a 


T 


Sec 0 
aO 


(30) 


where = zenith atmospheric transmission. Zenith atmospheric transmission 
ranges from 0.75 at 0.53 x 10"® meters to 0.83 x 10“® at 0.69 x 10"® meters 
under clear sky conditions but can of course drop to nearly zero under poor 
visibility conditions. A value of 0.70 is normally used in calculations. The path 
parameter can easily be evaluated on a desk ti^e computer. A solution for an 
orbital height of 1000 km is shown in Table 4. 
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Table 1 

Values of the Constant p for NASA 
Cube Corner Arrays 


Satellite 

P 

BE-B 

BE-C 

GEOS-I 

GEOS-n 

GEOS-C 

Timation III 

2.16 X lO'^ 

2.16 X 10“ 

2.00 X 10“* 

4.55 X lO'^ 

Not tested at writing date 
Data not reduced 


Table 2 

Maximum Velocity Aberration as a Function 
of Orbital Altitude 


Altitude (Km) 

Velocity Aberration (Rad x 10"®) 

0 

53 

500 

51.3 

1,000 

49.5 

2,000 

46.4 

4,000 

41.7 

8,000 

35.4 

16,000 

28.4 

32,000 

21.7 

64,000 

16.0 

128,000 

11.6 

256,000 

1 8.29 

512,000 

' 5.90 
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Table 3 

Variation of Apparent Velocity as a Function of Zenith Angle { 6 ) 
and Velocity Vector (w) for Several Altitudes 


s (Deg) 

<w(Deg) 

h = 

= 0 

h = 1000 Km 

h = 10,000 Km 
! 

Wv 

?7(Deg) 

V'/V 

^(Deg) 

V'/V 

T] (Deg) 

0 

0 

1.00 

0 

1.00 

0 


0 

15 

0 

1.00 

0 

1.00 

0 

1.00 

0 


45 

0.98 

44 

0.99 

44 

mmm 

45 


90 

0.97 

90 

0.97 

90 

rail 

90 

30 

0 

1.00 

0 


0 

1.00 

0 


45 

0,94 

41 


42 


45 


90 

0.87 

90 


90 

ilB 

90 

45 

0 

1.00 

0 


0 

1.00 

0 


45 

0.87 

35 


38 

mSm 

44 


90 

0,50 

90 


90 


90 

60 

0 

1.00 1 

0 


0 

1,00 

0 


45 

0.79 

27 

0.85 

34 

msm 

43 


90 

0.50 

90 

0.66 

90 


90 

73 

0 

1.00 

0 


0 

1.00 

0 


45 

0,73 j 


0.81 

29 

0.96 

43 


90 

1 

0.26 


0.33 

90 

0.92 

90 
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Table 4 

Satellite Cross Section at Zenith 



— 

Range 
(m X 106 ) 

Max Velocity 
Aberration 
(Bad X IQ-'G ) 

Effective 
Area* 
(cm^ ) 

■“1 

p Value 
(X lo'')* 

Gain 

(dB) 

Cross 
Section 
(m* X 10*)** 

Cross Section/ (Range) 4 

1. BE-B 

1.13 

49 

142 

2,16 

85,1 

4.60 

2.82 X 10"i® 

2. BE-C 

leOO 

49 

142 

2.16 

85.1 

4.60 

4.60 X 10'*® 

3. GEOS I (A) 

1.95 

46 

1793 

2,00 

85,0 

57.2 

3.96 X 10“*® 

4. GEOS n (B) 

1.53 

48 

2147 

4.55 1 

86.7 

100 

18.2 X 10'*® 

5. GEOS ni (C) 

0.93 

50 






6. TIMATION ni 

14.0 

29 

540 


92.8 

103 

1.62 X 10-2 

7. LAGEOS 

5.90 

38 

257 


86.0 

10.8 

8.91 X 10'^* 

8. LUNAK ARE AYS 

360 

7,0 

-^000 


96.0 

400 

2.4 X 10"2* 

9. STARLET 

0.92 

50 

101 


77.4 

0,55 

1 

7.67x10"*® 

1 


Effective Area and p value are computed for the satellite ot zenith. The p value is used to compute the gain function by the following relation 

G(si) =4vrP(v’;) p^e'^'f' =2p2 
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Figure 1 . 
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Figure 2a. Optical schematic of test setup for measuring cube array area. 
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Figure 2b. Optical schematic of test setup for measuring point spread function, 




AJAq - RELATIVE AREA 



Figure 3. Relative area of GEOS A cube corner panels as a function of incidence angle 


a- ARC SECONDS 


Figure 4. Far field patterns for cube corners used on several satellites*, 


ALL TESTS MADE AT NORMAL INCIDENCE 
ALL CUBE CORNERS HAD 2.54 CM HEXAGONAL 
ENTRANCE PUPILS 
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Figure 5. Coordinate system for far -field intensity pattern. 
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